Characteristics of Main Sequence Stars 


Main-sequence stars obey several relations (which are mostly pre- 
dictable from homology). 


e Main sequence stars obey a mass-luminosity relation, with £ œx 
M". The slope 7 changes slightly over the range of masses; be- 
tween 1 and 10Mo, 7 ~ 3.88. The relation flattens out at higher 
masses, due to the contribution of radiation pressure in the cen- 
tral core. (This helps support the star, and decreases the central 
temperature slightly.) The relation also flattens significantly at 
the very faint end of the luminosity function. This is due to the 
increasing important of convection for stellar structure. 


e Main sequence stars obey a mass-radius relation. However, the 
relation displays a significant break around 1Mo; R x M§, with 
€ = 0.57 for M > 1Mo, and € ~ 0.8 for M < 1Mo. This 
division marks the onset of a convective envelope. Convection 
tends to increase the flow of energy out of the star, which causes 
the star to contract slightly. As a result, stars with convective 
envelopes lie below the mass-radius relation for non-convective 
stars. (This contraction also increases the central temperature 
(via the virial theorum) and also moves the star above the nominal 
mass-luminosity relation. 


e Near the surface of a convective star, radiation often plays a 
role in the energy transport, so that the true value of the adia- 
batic gradient is Vraa > V > Vaa. Solving for V is tricky, as it 
involves a number of (poor) assumptions, and the adoption of a 
free parameter called the “mixing length”, which is the distance 
a blob of convective material moves before dissolving into its sur- 
roundings. (Actually, the mixing length is usually quoted as the 
ratio of this distance to the local pressure scale height (i.e., the 


length it takes the pressure to decline by 1/e), so that 


œ =lm/Ap where Ap=— 


Changing the mixing length changes V which changes the pressure 
gradient, which changes the size (and therefore the temperature) 
of the star. The larger the surface convective zone (or the smaller 
the pressure gradient), the larger the change in radius. 


e The mixing length is essentially a free parameter in the mod- 
els for stars with convective envelopes. For main sequence stars, 
uncertainties in the mixing scale length do affect the computed 
stellar radius (and effective temperature), but only slightly. A 
factor of two increase in a translates into a change of only +0.03 
in log Tig. The star’s luminosity, of course, is not affected. First 
generation stellar models (from the 1960’s) adopted a ~ 1, but 
those models produced stars that were too cool. A value of 1.5 or 
1.6 does a better job at fitting the observations. 


e The depth of the convective envelope (in terms of Meny/ Mr) 
increases with decreasing mass. Stars with M ~ 1Mo have ex- 
tremely thin convective envelopes, while stars with MS 0.3Mo 
are entirely convective. Nuclear burning ceases around M ~ 
0.08Mo. 


e The interiors of stars are extremely hot (T > 10° K). The fall-off 
to surface temperatures (T ~ 10* K) takes place in a very thin 
region near the surface. 


e The region of nuclear energy generation is restricted to a very 
small mass range near the center of the star. The rapid fall-off of 
€n with radius reflects the extreme sensitivity of energy generation 
to temperature. 


e Stars with masses below ~ 1Mo generate most of their en- 
ergy via the proton-proton chain. Stars with more mass than this 
create most of their energy via the CNO cycle. This changeover 
causes a shift in the homology relations for the stellar interior. 


e CNO burning exhibits an extreme temperature dependence. 
Consequently, those stars that are dominated by CNO fusion have 
very large values of £/4mr? in the core. This results in a large 
value of Vraa, and convective instability. In this region, convective 
energy transport is extremely efficient, and V ~% Vaa. 


e Because of the extreme temperature sensitivity of CNO burning, 
nuclear reactions in high mass stars are generally confined to a 
very small region, much smaller than the size of the convective 
core. 


e As the stellar mass increases, so does the size of the convective 
core (due again to the large increase in € with temperature). Su- 
permassive stars with M ~ 100Mo would be entirely convective. 


e A star’s position on the ZAMS depends on both its mass and its 
initial helium abundance. Stars with higher initial helium abun- 
dances have higher luminosities and effective temperatures. This 
is predicted by homology; the higher mean molecular weight trans- 
lates into lower core pressures. Helium rich stars therefore are 
more condensed, which (through the virial theorum) mean they 
have higher core temperatures and larger nuclear reaction rates. 


e Changes in metallicity shift the location of the ZAMS; the metal- 
poor main sequence is blueward of the solar metallicity main se- 
quence. This is primarily due to the reduced amount of bound-free 
absorption throughout the star (which only comes from metals). 
The smaller opacity allows the energy to escape more easily, and 


thus increases the luminosity via (3.1.4), i.e., 


16racr? T? dT 
3Kp dr 


L=- 


(The effect is slightly less for higher mass stars, since the amount 
of energy generation is proportional to Z in CNO burning, but the 
smaller opacity is still the most important factor.) The luminosity 
increase also results in a hotter surface temperature, via £L œx 
eT 


e Changes in the metallicity and helium abundance also affect the 
precise mass where the convective core disappears. Greater helium 
abundance forces pT to be larger to maintain the pressure. This 
produces an increase in luminosity, a larger value of Vyaq, and an 
increased tendency towards convection. Similarly, a decrease in Z 
translates into a lower opacity, a larger luminosity, and a greater 
probability of convection. 


e The central pressure is inversely proportional to stellar mass, 
as predicted by homology. Similarly, the central temperature also 
behaves as predicted, t.e., it is directly proportional to mass. At 
large and small masses, the central density also obeys its homology 
relation, but the central density becomes almost independent of 
mass near ~ 1Mo. This behavior is a natural consequence of the 
virial theorum, the main-sequence mass-luminosity relation, and 
the different temperature sensitivities of the CNO and pp chains. 


Since main sequence stars are in virial equilibrium, 


M M\° 
T, x — = T’ — eM? 18.1 
x o ( a x ppM (18.1) 
Meanwhile, the nuclear energy generation rate is 
dL 
— = En X poly = Lx p IY% M (18.2) 


dM 


while the mass-luminosity relation is 
L x M” (18.3) 


Now, if we substitute in for temperature using (18.1) and lumi- 
nosity using (18.3), then 


M" ox peTZ M 0 pe (3M2) M 


or 
3n—3—2v 


pe x M 3 (18.4) 


At high masses, v ~ 18 and 7) ~ 3, so the exponent is negative; at 
low masses; v ~ 5, and 7 ~ 2, and again, the exponent is negative. 
However, around 1Mo, v ~ 4, and 7 ~ 4; this generates a positive 
exponent, and thus changes the behavior of the central density. 
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The mass-luminosity relation for main-sequence stars. 
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The mass-radius relation for main-sequence stars. 
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The mass fraction of main sequence stars plotted against the total 
stellar mass. “Cloudy” areas indicate the extension of convection 
zones. The solid lines give the mass values at 1/4 and 1/2 the 
stellar radius. The dash lines show the mass within which 50% 


and 90% of the stellar luminosity is produced. 


The fraction of CNO energy generated for stars of different masses. 
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The central temperature and density for main sequence stars. 


Main Sequence Luminosity Evolution 


As a star burns on the main sequence, its luminosity will increase. 
This is due to the star’s higher mean molecular weight: since the 
core has fewer particles to support it, pT must increase, and so 
must the rate of nuclear reactions. To estimate the amount and 
timescale for this brightening, first consider that according to the 
virial theorum 


or 


T œ uM? pt (18.5) 


where we have substituted for R using p x M/R°. Now if we 
assume that energy is transported via radiation, then from (3.1.4), 


16acr? T? dT 


Pre 3Kp dr 


Very roughly, 


SO 


L x — (18.6) 


For main sequence stars, Kramers law opacity dominates (except 
perhaps in the very center, where electron scatter is important). 
Thus, 


RT15/2 


LX rs 


If we again substitute for R using p x M/R, and T using (18.5), 
then 
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(18.7) 
Since the mass of the star is (essentially) constant, and the density 
dependence is very weak, (18.7) says that 


(18.8) 


where w = 15/2. (Note that the homology equation for a Kramer 
law opacity gives gives 8.2 < w < 9.0, depending on the pre- 
cise value of v.) Next, consider the definition of mean molecular 
weight. From (5.1.5) and (5.1.6), the mean molecular weight of 


the electrons is 
Tuy 2 
n= (SBA) =e 


Meanwhile, the mean molecular weight of ions is 
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which reduces to 
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The total mean molecular weight, as a function of hydrogen frac- 
tion, is therefore 


HI (18.9) 


MOR | ee | : (18.10) 


Finally, if the fusion process releases Q = 6 x 1018 ergs for every 
gram of hydrogen fused, then the rate of fuel consumed is related 
to the star’s luminosity by 


255, ee (18.11) 


Through equations (18.8), (18.10), and (18.11), we can compute 
the how luminosity of a main sequence star changes with time, 
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Putting it all together, the expression for £ becomes 
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This is easily integrated to give 
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If uo ~ 0.6 and Y% = 15/2, then in solar units, we have 
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For the sun, this relation gives Lo œ~ 1.26L 0; if the homology 
exponents are used, the value changes to between 1.28 and 1.31 
This compares to detailed models which imply that the sun has 
brightened by a factor of 1.37. (In other words, the sun asa ZAMS 
star was ~ 25% less luminous than it is today.) 


Other Facets of Main Sequence Evolution 


e During main sequence evolution, the increase in core tempera- 
ture increases the pressure just outside the core. This (through 
some complicated physics), causes the surrounding regions to ex- 
pand. (This in part explains why our estimate of the solar lu- 
minosity increase overestimates the actual number. Some of the 
energy has gone into expansion.) In general, any time a stellar 
core contracts, the envelope surrounding the core will expand. 


e In high mass stars, the CNO nuclear reactions are highly concen- 
trated in the center of the star. As a result, the molecular weight is 
unchanged over most of the star, and only the core contracts. The 
star therefore becomes redder, as the increase in the radius over- 
whelms that of the luminosity. For low mass stars, the relatively 
weak dependence of pp-cycle energy generation on temperature 
allows more of the star to be involved in nuclear reactions. A 
larger fraction of the star thus has its molecular weight changed, 
and more of the star participates in the contraction. This com- 
petes with the overall stellar expansion. (A 1M o main-sequence 
star actually becomes hotter with age.) 


e As a fully convective star (M < 0.3Mo) evolves, the change 
in its molecular weight just moves the star to a higher helium 
content ZAMS model. 


e Convective cores in stars with M < 2.25Mo are small, but 
they increase (slightly) with time. As the core temperature in- 
creases, the CN-cycle becomes more important, driving more con- 
vection. Eventually, however, the convective core will again be- 
come smaller, as the decrease in Kramer’s opacity takes over. 


Semi-Convection 


For stars with M > 2.25Mo, the mass fraction of the convective 
core decreases with time. ‘This is due to the decrease in the opac- 
ity with temperature, which follows from Kramer’s law. As the 
convective core gets smaller, it leaves behind a chemical gradient. 
This causes a phenomenon known as semi-convection. 


Consider a star with a convective core within which nuclear reac- 
tions are proceeding, and a radiative envelope. In the convective 
core, Vraa > Vaq; outside the core, Vraa < Vaq, and, at the 
border between the two regions, Vraa = Vaa. AS always, Vraq is 
given by 


3kKLP 
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Now assume that the major source of opacity is electron scat- 
tering, and that the matter is fully ionized. (Since we are talking 
about the core of a massive star, these are excellent assumptions. ) 
As nuclear reactions proceed in the core, the hydrogen content of 
the region will decrease. Since the electron scattering opacity for 
fully ionized matter is given by 


Mradi = (3.1.6) 


Ke = 0.2(14+ X) (6.1.7) 
this results in a decrease in the stellar opacity. 


Now let’s examine the behavior of V;aq at a distance € from the 
convective core. By definition, for convection to occur 


Vradi > Waa inside the core 


and 
Vrado < Vaq outside the core 


Therefore, since Vaq shouldn’t change much over a distance e€, 
Vradi > Vrad,o by the definition of convection. However, Vrad 
is directly proportional to the hydrogen content, and X; < Xo. 
Thus, Vrad,i < Vradjo from the opacity. In other words, there is a 
fundamental contradiction. 


Another way to look at the problem is to consider the temper- 
ature and density difference between a convective blob and its 
surroundings. Recall that for a region to be dynamically stable, 


DO. 


or 


where A represents the difference between the region internal to 
the blob and the blob’s surroundings. For positive displacements 
of Ar, the blob will be denser than its surroundings, and will sink 
back to its original position. 


Now consider the temperature difference between the blob and 
the surrounding. 


s- (2) 


If the blob is in pressure equilibrium (given) and the stellar region 
is chemically homogeneous, then from (3.2.2) 


The negative sign indicates that the inside of the blob will be 
cooler than its surroundings. It will therefore gain some energy 


by radiation, which will result in the blob having a slightly lower 
density, and a slightly weaker restoring force due to buoyancy. 
Any motion of the blob will eventually be damped out. 


Next, consider the situation where the surrounding medium is 
chemically inhomogeneous. We again assume that the region is 
dynamically stable, with Ap/Ar > 0. However, in this case 


= +E (8.2.1) 


Consequently, if 


d d A 
Ay = (H) — (4) > dex 
dr); dr } . p 


then blobs that have a positive Ar displacement will have AT > 
0. These blobs will be hotter than their surroundings and will 
therefore lose heat as they move. When this happens, their density 
will increase, the restoring force will increase, and their oscillation 
will slowly increase in amplitude. The blob will be vibrationally 
unstable; the result will be a slow mixing, which will eventually 
destroy the chemical gradient. This semi-convection will occur 
when 


p (dlou 
Vad < V < Vaat 5 (T) (8:2:2) 


Note that from (8.2.1) and (8.2.2), it will only occur when the 
mean molecular weight is decreasing rapidly with radius, or equiv- 


alently, 
dln u 
€ In = | i 


This occurs in high mass stars with convective cores. In these ob- 
jects, the convective region gets smaller with time, leaving behind 
a chemical gradient. 


